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Abstract. 

■ The entanglement of two dipole-coupled atoms with nondegenerate two-photon 

transitions interacting with two-mode field in lossless cavity has been investigated. It 
shows that the entanglement is dependent on the initial atomic states. The possibility 
of considerable growth of atomic entanglement due atomic coherence and dipole-dipole 
^ ' interaction is shown in the case of great mean values of thermal photons. 
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1. Introduction 

Entanglement is a key resource whicli distinguislies quantum from classical information 
theory. It plays a central role in quantum information, quantum computation and 
communication, and quantum cryptography [l]. In order to function optimally these 
applications require maximally entangled states. However, interactions with the 
environment always occur, and will degrade the quality of the entanglement. Although 
the interaction between the environment and quantum systems can lead to decoherence, 
it can also induce entanglement [2] . Thus, understanding and investigating entanglement 
of mixed states becomes one of the actual problem of quantum information. Recently, 
Bose et al. |3j have shown that entanglement can always arise in the interaction of 
an arbitrary large system in any mixed state with a single qubit in a pure state, and 
illustrated this using the Jaynes-Cummings interaction of a two-level atom in a pure 
state with a field in a thermal state at an arbitrary high temperature. Kim et al. |1] 
have investigated the atom-atom entanglement in the system of two identical two-level 
atoms with one-photon transition induced by a single-mode thermal field. They showed 
that a chaotic field with minimal information can entangled atoms which were prepared 
initially in a separable state. Zhou et al. [6] have considered the same problem for 
nonidentical atoms with different couplings. The entanglement between two identical 
two-level atoms through nonlinear two-photon interaction with one-mode thermal field 
has been studied by Zhou et al. [5]. They showed that atom-atom entanglement 
induced by nonlinear interaction is larger than that induced by linear interaction. In [7] 
has discovered that two atoms can be entangled also through nonlinear nondegenerate 
two-photon interaction with two-mode thermal field. The influence of dipole-dipole 
interaction on entanglement between two cubits induced by one-mode and two-mode 
thermal field has been investigated in [8]-|10]. 

The problem of creating or controlling the atomic entanglement is greatly related 
to the atomic coherence of population between different levels. Hu et al |II],[12] have 
shown that the entanglement between two atoms induced by one-mode thermal field 
can be manipulated by changing the initial parameters of the atoms, such as the 
superposition coefficients and the relative phases of the initial atomic coherent state 
and the mean photon number of the cavity field. In [13j Hu and Fang study the 
effect of the atomic coherence on the entanglement of two two-level atoms interacting 
with two-mode thermal fields through a nondegenerate two-photon process, which 
is a nonlinear interaction They have found that for some atomic initial state the 
entanglement induced by nonlinear interaction may be larger than that induced by linear 
interaction. They have also discovered that entanglement may be greatly enhanced due 
to atomic coherence. This provides a method by which one can control the entanglement 
of the whole system. But the influence of the dipole-dipole interaction on the coherence- 
enhanced entanglement has not be considered. 

In this paper we study the effect of the atomic coherence on two atoms entanglement 
induced by two-mode thermal fields taking into account the dipole-dipole interaction. 
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In this model, we consider a system composed of two identical two-level atoms coupled 
via dipole-dipole interaction. The frequency of the atomic transition for atoms is uiq. 
Both atoms are also interacting with two modes of quantum electromagnetic field with 
frequencies uii and ijJ2 through a nondegenerate two-photon process. For simplicity 
we ignore the Stark shift and assume the lossless cavity modes are at the two-photon 
resonance with the atomic transition, i.e. the condition uji+liJ2 — ojq takes place. Then, 
the Hamiltonian of the considered system in interaction picture can be written in the 
following form: 



(1) 



i=l 



where and aj are the creation and the annihilation operators of photons of jth cavity 
mode {j — 1,2), and R^ are the raising and the lowering operators for the ith 
atom {i — 1,2), g is the coupling constant between atom and the cavity field, Q is the 
coupling constant of the dipole interaction between the atoms and |-|-) and |— ) are the 
excited and the ground states of a single two-level atom. The two-atom wave function 
can be expressed as a combination of state vectors of the form \ vi,V2) — \ vi\\v2), where 

VuV2 = +, -. 

The density operator for the atom-field system follows a unitary time evolution 
generated by the evolution operator Ui{t) = exp[—tHit/h]. On the two-atom basis, 
1+)+)) analytical form of the evolution operator Ui{t) is 

given by 



U{t) 



( Un Uu C/i3 Uu \ 

U21 U22 U23 U24 

U31 U32 U34 

U41 U42 ^43 U44 



where 



AAA 

Uii — 1 + 2aia2 Y<^i^<^2^ ) Uu = 2aia2-raia2, U41 = 2afa2Y^t(^tj 
AAA 



A 



B 



B 



C/44 = 1 + 2a'{a2—aia2, U12 = t/13 = aia2—, U21 = t/31 = -j^i ai , 

B B 

U24 = U34 = —0,10,2, U42 = U4S = 0^02 



U22 — 



exp 



-i^{a + e)t 



(2) 



|[1 - exp{ig9t)]a + 2^exp(^|(3a + 9)t] + 9[1 + exp(«^^t)]| , 
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C/23 = u. 



32 



exp 



-i^{a + e)t 



AO 



[1 - e-x.p{iget)]a - 2^exp(z|(3Q; + e)t\ + e[l + exp(z(/^i)] [> , 



and 



A = exp 



9(^. 
-I — t 
2 . 



B — exp 



2 



[1 - exp{igdt)] , 



a — —, X — 2(aia2afa2 + aia2aia2), 6 — \/8(aia2aia2 + aia2aia2) + a"^- 
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The initial cavity mode state are assumed to be the thermal two-mode state 



Pf(0) = X^H^'i (^1)^2(^2) |ni,n2)(ni,n2|. 
The weight functions are 



(3) 



m n2 



Pi{ni) 



1 + 



\ni+l ■ 



where fii is the mean photon number in the ith cavity mode, fti = {exp[huJi/kBT] — 
kB is the Boltzmann constant and T is the equilibrium cavity temperature. 

We consider also the initial state of each atoms to be prepared in a coherent 
superposition of the two levels, that is, 

l^i(O)) = cos^i|+) + e*^isin^i|-), |^2(0)) = cos^aK) + e*^^ sin^sh)- 

Here 61 and 62 denote the amplitudes of the polarized atoms, and and ip2 are relative 
phases of the two atoms, respectively. So the initial density matrix for the two atoms 
can be written as 



Pa{0) 



V 



/ 



(4) 



/ Pll(O) Pl2(0) Pl3(0) pi4(0) \ 

PuiO) P22(0) P23(0) P24(0) 

PUQ) PUQ) P33(0) P34(0) 

Pu(0) ^^4(0) P^4(0) P44(0) 

where the matrix elements are expressed as follows: 

Pii(O) = cos^ 9i cos^ 02, Pi2(0) — cos^ 6*1 cos ^2 sin^2e~'^^, 
Pi3(0) = cos^isin^iCOs2^2e~"^\ pu{0) = cos ^1 sin ^1 cos ^2 sin ^2e~'^'^^+'^'\ 
P22(0) = cos^ 01 sin^ 02, P23(0) = cos 01 sin 9i cos 02 sin ^26"'^'^^"'^^^ 
P33(0) = sin^ 01 cos^ 02, P24(0) = cos 01 sin Oi sin^ ^26""^^ , P34(0) = sin^ 61 cos 02 sin 02e-"^\ 

P44(0) = 1 - Pll(O) - P22(0) - P33(0). 
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3. Results 

To investigate the entanglement between atoms one can obtain the time- dependent 
reduced atomic density operator by tracing the combined atom-field density operator 
over the field variables: 

PAit) = TrpU{t)pF{Q) ® pA(0)f/+(t). (5) 

For the two-qubit system described by density operator pA{t), a measure of entanglement 
can be defined in terms of the negative eigenvalues p^ of partial transpose [ll],[l5] of 
the reduced density matrix 

e = -2j2p-. 

When 5 = two atoms are separable and e > means the atom-atom entanglement. 
The case e = 1 indicates maximum entanglement. 

Substituting Eqs (2)- (4) into Eq. (5) one can obtain the atomic density matrix at 
time t, accordingly, we can write down the partial transpose matrix 



PA{t) 



( Pll 



Pl2 Pl3 

P22 P23 

P23 P33 

P24 P34 



Pl4 \ 
P24 
P34 
P44 y 



(6) 



the evident expressions of the matrix elements are given in Appendix. 

In Fig. 1 we plot the negativity as a function of gt for a fixed value of dipole 
strength a = 0.1 and weak thermal two-mode field with ni = 77-2 = 0.01. The Fig 
la) corresponds to atoms in the incoherent states and Fig.2b) corresponds to them in 
two different coherent states. Comparing Fig. 1 with that for atoms without dipole 
interaction [13], we find that they are clearly different. Ignoring the dipole interaction 
we can obtain that the atom-atom entanglement is greatly enhanced owing to the atomic 
coherence. On the contrary the presence of an atomic coherence for dipole coupled atoms 
leads to reduction of the entanglement maximum value. In Fig. 2 we plot the negativity 
as a function of gt for a fixed value of dipole strength a = 0.1 and weak thermal two- 
mode field with ni = = 0.2 but different values of relative phase A(f = (fi — ip2- The 
atoms are assumed to be prepared in coherent states with ^1 = 6*2 = 7r/4. Clearly, the 
entanglement may be stronger for appropriate value of relative phase Ayj. In Fig. 3 and 
Fig. 4 we show the entanglement time dependent for intensive two-mode thermal field. 
With increase of the mean photon numbers (see Fig. 1 - Fig. 4) the value of atom-atom 
negativity decreases. But comparing Fig. 3 and Fig. 4 with the same in [13] one can 
easily obtain that in view of dipole interaction we gain essentially large enhancement 
of entanglement degree transferring from incoherent to coherent atomic initial state. 
The results indicate that we can produce more strong entanglement by initial atomic 
coherence and dipole-dipole interaction. 
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a) b) 




Figure 1. The negativity as a function of gt for the model with m = n2 = 0.01 and 
a = 0.1: a) = |+), |*(0))2 - |-), b) = (1/72(1+) + h)), |*(0))2 = 

(1/V2(|+) - I-)) (soUd) and |vl/(0))i = (1/^2(1+) + |-)), |*(0))2 = (l/v^(|+) + |-)) 
(dashed). 

a) b) 



e(t) e{t) 




Figure 2. The negativity as a function of gt for the model with m = n2 = 0.2 
and a = 0.1, and atoms to be prepared in coherent states |^'(0))i — (l/-\/2(|+) + 
|-)),|^-(0))2 = (1/V2(|+) + I-)): a) A(p = 0, b) A<^ = tt (solid) and Atp = 7r/6 
(dashed). 



4. Conclusion 

In this paper we have investigated the effect of the atomic coherence on the entanglement 
of two two-level identical atoms interacting with two-mode thermal field through a 

nondegenerate two-photon processes in the lossless cavity taking into account the dipole 
interaction. It is shown that the entanglement is dependent on the initial atomic 
state and strength of the dipole interaction. The results also show that the atom- 
atom entanglement can be controlled by changing the system parameters, such as 
relative phases, the amplitudes of the polarized atoms, the mean photon numbers of 
the two-mode thermal field, and the strength of the dipole interaction (or distance 
between atoms). The have derived that the dipole-dipole interaction can produce the 
appreciable amount of entanglement even in the case of large thermal noise if the atoms 
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a) b) 



e(t) e{t) 




Figure 3. The negativity as a function of gt for the model with a = 0.1: a) ni = 

n2 = 10, |vl/(0))i = |+),|*(0))2 = I-) (dashed), \<\>{Q)), = (1/V2(|+) + |*(0))2 = 
(l/^/2(|+) - I-)) (soUd), b) ni = n2 = 40, |*(0))i = |+), |^'(0))2 = |-) (dashed), 
If (0))i = (1/V2(|+) + I-)), 1* (0))2 = (1/V2(|+) - I-)) (sohd). 

a) b) 




Figure 4. The negativity as a function of gt for the model with m = n2 = 20: a) 
a = 0.3, |vl>(0))i = |+), |*(0))2 = I-) (dashed), |vl/(0))i = (1/^2(1+) + H), |*(0))2 = 
(1/^2(1+) - I-)) (solid), b) a = 1, |v|/(0))i = |+), |vE'(0))2 = I-) (dashed), |vl/(0))i = 
(1/72(1+) + I-)), |v&(0))2 = (1/V2(|+) - I-)) (solid). 



to be prepared in coherent states. In this paper we consider only one type of the 
atomic initial coherent states. The entanglement behavior for initial atomic coherent 
states |*(0)) = cos 611+, -)+ sin 611-,+) and |*(0)) = cos 6*1+, +) + sin 611-, -) will be 
considered in our following paper. 
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Appendix: Evident expressions for density matrix Elements in Eq. (6) 



oo oo / / A \ 

Pn{t) = E E Pi{ni)P2{n2)pn{0) 1 + 2{n, + l)(n2 + l)^^l±l^ 

ni=0n2=0 \\ ■^ni+l,n2+l / 



X 



A* \ I R 



2' 



X 1 + 2(ni + l)(n2 + l)^ii±i£ii±i + (^^^(O) + ^32(0) + ^23(0) + P33(0))mn2^i^ + 

\ A„i+i,n2+l/ t:'ni,n2 / 



l^ni — l,n2 — 1 1 



+4 E E Pi("'i)P2(?T'2)p44(0)nin2(ni - l){n2 - 1) 

711=2 n2=2 ^^ni-l,n2-l 

00 00 

P12(^) = E E Pl(m)p2(n2)(pi2(0)(C/2*2) ni,n2 + Pl3(0)(C/22)ni,n2) ^ 
ni=0 n2=0 

X f 1 + 2(ni + l)(n2 + + 

\ A„^+l,n2+l / 

00 00 TD TD* 

-Dni,n2 -'-'ni — l,n2— 1 



+ E E Pi(^i)P2(n2)(p24(0) + P34(0))nin2 



/3 /3 ' 

"ni,n2 ^ni—l,n2 — l 



ni = l 712 = 1 ^711,712 ^711-1,712- 

Pl3it)= E EPl(^l)P2(n2)(pi2(0)(t/3*2) 711,712 + Pl3(0)(t/33)ni,7l2) ^ 
711=0 712=0 

x(l + 2(ni + l)(n2 + + 



V ^ni+l,n2+l 

00 00 

+ E E Pl(^l)P2(w2)(P24(0) + P34(0))nin2 



R R* 

-0711,712 -^ni — 1,712-1 



711=1712 = 1 ^711,712 ^711-1,712-1 



00 00 / / ^ 



PuW -EE PiK)P2(n2)pi4 ( ( 1 + 2(ni + l)(n2 + lY-^ll±l^ X 

711=0 712=0 VV A„i + 1,„2 + 1 / 

X fl + 2{n, + l)(n2 + i)^^lzlmzl] + p,{0)p,{0)pu{0) (l + 2^]] + 

V /^Tll -1,742 -1 / \ ^1,1// 

,^7H + 1,1 \ I N / 1 I O ^1,712 + 1 



+ E Pl(m)p2(0)pi4 1 + 2-^^ + E Pl(0)p2(n2)pi4 1 + 2- 



711=1 V A„j+i,iy ^^^^ \ Ai+i,n2+l, 



P22(t) = E E Pi(ni)P2(n2) (pn(0)(ni + l)(n2 + 1) '^7^''"^''" + 

711=0 712=0 \ ^711 + 1,712 + 1 

+P22(0)(C/22)7^i,7.2(t^22):„„2 + P32(0)(C/23)7^l,7.2(t^22):„„2 + P23(0) (C/22)7.l,7^2 (t^23):„7.2 + 

I R |2 

|-'-'7ll — 1,712 — 1| 



|2 



00 00 



+P33(0)(t/23) 711,712 "1,712) + E E Pi (^1)^2(^2)^44(0)^1712- 2 

711=1 712 = 1 "711-1,712-1 



00 00 / I R |2 

P23{t) = E E Pi(^i)P2(n2) Pii(0)(ni + l)(n2 + l)^l±i^!^+ 



57tl+l, 712 + 1 I 

^2 

711=0 712=0 \ ^7il+l,7l2 + l 

+P22(0)(^722)77i,772(^32);i,n2 + P32 (0) (^733)771,772 (^32):i,„2 + P23 (0) (^722)711 ,„2 (t^33):i,772 + 
+P33(0)(C/23)77i,7.2(C/33):i,„2)+ E E Pl(^l)P2(n2)nin2P44(0)^, 

711=0 712=0 ^711,772 
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oo oo 



ni=0n2=0 ^ni,n2 ^ni+l,n2+l 



oo 



+ E Pl(^l)P2(0)(p24(0)(t/22)ni,0 + P34(0)(t/23)ni,0) + 



ni=l 



+ E Pl(0)P2(n2)(p24(0)(t/22)o,n2 +P34(0)(t/23)o,nJ + 
n2=l 

+Pl(0)p2(0)(p24(0)(C/22)o,0 + P34(0)(C/23)o,o) + 

oo oo 

+ E E Pi (^1)^2(^2) (P24(0) (C/22)„i,n2 +P34(0)(C/23)ni,n2)x 

ni=l n2=l 

X fl + 2nin2-^"^~''"^~'^ 



A, 



m— 1,712—1 



P33W E E Piini)p2MpniO) ( (n, + l)(n2 + i)\^!^l±h!^+ 

ni=0n2=0 \ 
+P23(0)(t/32)ni,n2(t^3*3)ni,n2 + P32(0) (C/33)ni,n2 (t^3*2)ni,n2 + 
+P22(0)(C/32) ni,n2 ni,n2 + P33(0)(C/33) ni,n2 ni,Ti2 ) + 

+ E E Pl(^l)P2(n2)p44(0)nin2 '''-^-"-'-^ ^ 
ni=ln2=l fni-l,n2-l 

00 00 JD* 75 

P34W E E Pi(m)P2(n2)pn(0)K + 1)(^2 + i)^^l±i£-±i+ 

ni=0n2=0 ^ni,n2 f'm+l,n2+l 

00 

+ E Pl(^l)p2(0)(p24(C^32)ni,0 + P34(C^33)ni,0+ 

ni=l 

00 

+ E Pl(0)p2(n2)(p24(C^32)o,n2 +P34(C^33)o,n2 + 

n2=l 

+pl(0)p2(0)(p24(t/32)o,0 + ^34(^/33)0,0+ 

00 00 

+ E E Pl(^l)P2(n2)(p24(0)([/'32)ni,n2 +P34(0)(t/33)ni,n2)X 
ni=l n2=l 

X fl + 2nin2-'^"^"''"^"'^ 



A, 



ni — l,n2— 1 / 
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